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SUMMARY The Drosophila trachealess gene encodes a 
basic-helix-loop-helix-PAS transcription factor that controls 
the formation of the trachea and salivary duct. An ortholog of 
trachealess was identified in the brine shrimp, Artemia fran- 
ciscana, and was shown to be highly conserved by sequence 
identity. Expression of Artemia trachealess was observed at 


two sites during development: the naupliar salt gland and the 
juvenile thoracic epipod. These two organs function at their 
respective times of development in osmoregulation, an im¬ 
portant aspect of brine shrimp physiology. This extends the 
range of putative functions of trachealess to include forma¬ 
tion of osmoregulatory, respiratory, and ductile organs. 


INTRODUCTION 

Osmoregulation is a fundamental aspect of organismal phys¬ 
iology. It is particularly important for branchiopod brine 
shrimp, such as Artemia franciscana. Artemia survive in the 
hypersaline environments of salt ponds and brine pits. A re¬ 
quirement for survival of halophilic animals is possession of 
efficient mechanisms for both osmotic and ionic regulation. 
Water uptake occurs in the alimentary tract by swallowing 
salt water, and the excess NaCl is actively excreted (Croghan 
1958b; Hootman and Conte 1974). Artemia deal with the 
problem of salt excretion using different organs at different 
stages of their life cycle. Nauplii, which are the early stages 
of larval development, rely on a transient structure, the salt 
gland, for salt excretion (Conte 1984), whereas juveniles and 
adults rely on the thoracic epipods (Croghan 1958a). The an¬ 
tennal glands and maxillary glands may also contribute to 
osmoregulation (Martin 1992). Artemia cysts possess a fully 
functional salt gland when they hatch into nauplii (Croghan 
1958a). Beginning at the juvenile stage, the salt gland is re- 
sorbed concomitant with the appearance of the thoracic epi¬ 
pods (Conte 1984). The cells of the salt gland and the epi¬ 
pods share morphological features consistent with their roles 
in electrolyte transport, including extensive membrane am¬ 
plification and a rich supply of mitochondria (Conte 1984). 
Relatively little is known about the molecular processes that 
underlie the formation, differentiation, and function of the 
Artemia osmoregulatory organs. However, the Artemia 
APH-1 gene is expressed in the salt gland (Chavez et al. 
1999). APH-1 encodes a POU-III subclass homeobox-con- 
taining transcription factor highly related to the Drosophila 


Drifter (Dfr) protein (Anderson et al. 1995). One prominent 
role of Drosophila Dfr is in controlling tracheal development 
and transcription by acting as a coactivator of the Trachea- 
less (D-Trh) transcription factor (Zelzer and Shilo 2000). 
Here we identify an Artemia trachealess ( Af-trh ) ortholog. 

The D-trh gene is expressed in the salivary tissue, trachea, 
posterior spiracle, and central nervous system (CNS) (Isaac 
and Andrew 1996; Wilk et al. 1996). D-trh is expressed in 
the salivary primordium, which includes the precursors to 
both salivary duct and glands, but later expression is re¬ 
stricted to only the salivary duct. Genetic analysis of D-trh 
mutant embryos revealed that D-trh is required for the for¬ 
mation of the salivary duct but not the salivary gland (Isaac 
and Andrew 1996). During tracheal development, D-trh is 
expressed in the tracheal placodes and continues to be ex¬ 
pressed in all tracheal cells. Mutations in D-trh result in an 
absence of tracheal cells and a reduction in expression of nu¬ 
merous genes that are normally expressed in the trachea 
(Younossi-Hartenstein and Hartenstein 1993; Isaac and An¬ 
drew 1996; Wilk et al. 1996). The role of D-trh as a master 
regulator of tracheal development is reinforced by the obser¬ 
vation that misexpression of trh can induce ectopic forma¬ 
tion of tracheal cells (Wilk et al. 1996). Analysis of trh in an¬ 
other insect, the silkworm Bombyx mori , shows that it is also 
expressed in the trachea, supracolonic trachea (analogous to 
Drosophila posterior spiracle), and the anterior region of the 
developing silk gland (Matsunami et al. 1999). The silk 
gland is derived from the labial segment, as is the Dipteran 
salivary primordium. Although Bombyx trh ( Bm-trh ) is ini¬ 
tially expressed in the entire elongating silk gland, expres¬ 
sion eventually becomes restricted to just the anterior silk 
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gland, which functions as a duct for silk protein transport to 
the spinneret. This is analogous to the restriction of expres¬ 
sion and function of Drosophila trh in salivary duct forma¬ 
tion (Isaac and Andrew 1996). 

The trh gene encodes a protein of the basic-helix-loop- 
helix-PAS (bHLH-PAS) family of transcription factors. This 
class of proteins is highly conserved throughout evolution 
and is involved in a variety of developmental and physiolog¬ 
ical processes in organisms as diverse as Caenorhabditis el- 
egans, Drosophila , and human (Crews 1998; Gu et al. 2000). 
bHLH-PAS proteins are characterized by a DNA binding 
basic region, HLH domain for dimerization to bHLH-PAS 
protein partners, and the PAS-1 and PAS-2 domains that are 
multifunctional interaction domains. The D-Trh protein func¬ 
tions as a heterodimer in vivo with the Tango (Tgo) bHLH- 
PAS protein (Ohshiro and Saigo 1997; Sonnenfeld et al, 
1997). Together, they bind to DNA sequence elements on tar¬ 
get genes and activate transcription. The sequence elements 
contain an ACGTG core sequence and are referred to as a CNS 
Midline Element (Wharton et al. 1994), because it is also a tar¬ 
get sequence of Single-minded (Sim):Tgo heterodimers, 
which control CNS midline cell development and transcrip¬ 
tion. Activation of tracheal gene expression by D-Trh;Tgo re¬ 
quires a direct interaction between D-Trh and Dfr via the 
D-Trh PAS domain (Zelzer and Shilo 2000). 

Given the prominent roles that D-trh plays in development 
of multiple Drosophila tissues and the existence of a D-Trh co¬ 
activator in Artemia (APH-1), we have asked whether there is 
an Artemia trh gene. The results demonstrate the presence of an 
Artemia gene ( Af-trh ) highly related to insect trh that is ex¬ 
pressed in the salt gland and thoracic epipods. Although mor¬ 
phologically and developmentally distinct, the salt gland and 
epipods both function in osmoregulation, suggesting that Af-trh 
regulates genes involved in that process. This expands the po¬ 
tential developmental roles of trh to include osmoregulatory or¬ 
gans, in addition to controlling insect tracheal, salivary duct, 
and silk gland development. Because APH-1 is also expressed 
in the salt gland, this suggests that Af-Trh and APH-1 may in¬ 
teract to control salt gland transcription, similar to their interac¬ 
tion that controls tracheal transcription in Drosophila. 

MATERIALS AND METHODS 
Artemia culture 

Artemia franciscana were purchased from Carolina Biological. 
Cysts (Burlington, NC, USA) were hydrated and cultured in recon¬ 
stituted seawater at 20°C. 

Construction of Artemia genomic library 

Artemia trh genomic clones were isolated from a library of bacte¬ 
riophage X clones containing Artemia genomic DNA. Artemia were 
grown for 48 h and collected using a Nitex filter (Tetko, Lancaster, 


NY, USA). DNA isolation was performed using the Easy-DNA kit 
(Invitrogen, Carlsbad, CA, USA) and subjected to partial digestion 
with Sau3A. Sucrose gradient fractionation of the resulting frag¬ 
ments was carried out, and fractions containing DNA 20-25 kb in 
size were collected and pooled. These DNA fragments were ligated 
into the vector XDashll/SamHI. The recombinant clones were 
packaged, plated, and amplified. 

Isolation and analysis of Af-trh genomic clones 

Af-trh genomic clones were identified using a 3: P-labeled probe con¬ 
taining the bHLH region of the D-trh gene (Isaac and Andrew 1996; 
Wilk et al. 1996). The probe was generated by polymerase chain re¬ 
action (PCR) amplification of a D-trh cDNA clone, provided by 
B.-Z. Shilo (Weizmann Institute of Science, Rehovot, Israel), using 
the primers, 5'-CGCAAGGAGAAGTCGAGG-3' and 5'-CAGCT- 
TCAGGTAQGCTGAT-3', corresponding to amino acids 78-129. 
Hybridization to fdters containing 1.5 X 10 ■ recombinant phage 
was carried out in 50% formamide, 6X SSPE, 0.1% NaPPi, 0.1% 
SDS, 5X Denhardf s, 100 fig/ml salmon sperm DNA at 42°C for 18 h. 
After hybridization, filters were washed three times for 15 min each 
in 2X SSC, 0.1% SDS at room temperature and two times for 30 
min each in 2X SSC. 0.1% SDS at 50°C. DNA from positive clones 
was isolated and subjected to restriction enzyme cleavage, followed 
by Southern blot analysis and hybridization to the D-trh bHLH 
probe. Af-trh fragments were subcloned and sequenced. 

Generation and analysis of an Af-trh cDNA clone 

Reverse transcriptase (RT)-PCR was used to generate an Af-trh 
cDNA clone using RNA from animals collected 1 to 14 days after 
hydration, Poly(A)+ RNA was purified using a Poly(A) Pure kit (Am- 
bion, Austin, TX, USA). cDNA synthesis was performed using a 
primer containing oligo(dT) and additional 5' sequences for subse¬ 
quent PCR (5'-CCAGTGAGCAGAGTGACGAGGACTCGAGC 
TCAAGC[T18]-3') with 5 ng of Poly(A)+ RNA. PCR was carried 
out on cDNA in two steps as a nested reaction. Primers for the 
first PCR were (1) 5' -T ATTTTTTC AGT ATTTT GG A A ATG -3' 
and (2) 5'-CCAGTGAGCAGAGTGACG-3\ Primers used in the 
second PCR were (1) 5' -GATC AAGGCGTGG A A A AG-3' and (2) 
5' -GAGGACTCGAGCTCAAGC-3'. The PCR amplification prod¬ 
ucts were subcloned into pCRH-Topo-TA (Invitrogen) and sequenced. 

Sequence analysis 

The conceptual Af-Trh protein sequence was aligned to existing 
bHLH-PAS proteins using CLUSTAL W (Thompson et al. 1994). 
The limits of the PAS-1 and PAS-2 domains correspond to regions 
defined by Taylor and Zhulin (1999). Phylogenetic trees were gen¬ 
erated from CLUSTAL W alignments of bHLH-PAS proteins using 
the neighbor-joining algorithm (Saitou and Nei 1987). 

Whole mount in situ hybridization 

The in situ hybridization protocol was modified from Nulsen and 
Nagy (1999). Animals were collected from 1 to 14 days after hydra¬ 
tion and fixed in 4% paraformaldehyde, 1X phosphate-buffered sa¬ 
line for 2 h. They were then processed through an ascending meth¬ 
anol series (25%, 50%, 75%, 100%) and stored at 4°C after washing 
four times in 100% methanol. Preparation for hybridization in¬ 
volved rehydrating animals through a decreasing series of methanol 
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into PBT (1X phosphate-buffered saline + 0.1% Tween). After 
washing four times in PBT, animals were sonicated (Branson probe 
tip sonicator setting at 1.5, Danbury, CT, USA) until approximately 
50% were visibly destroyed (1-2 min). The sonication step en¬ 
hances probe accessibility. Artemia were then washed once in PBT 
and three times in 0.1 m triethanolamine. Even in the best prepara¬ 
tions, only 50% of the animals showed positive staining. 

Digoxygenin (DIG)-labeled RNA antisense and sense probes 
were synthesized corresponding to a 1.5-kb Af-trh cDNA fragment 
and 1.2-kb Artemia engrailed ( Af-en ) cDNA fragment. These probes 
were synthesized using either T7 or T3 RNA polymerase (Promega, 
Madison, WI, USA) from a linearized pBluescript II template con¬ 
taining either an Af-trh or Af-en cDNA clone fragment. Animals 
were incubated at room temperature for 20 min in triethanolamine 
with acetic anhydride (1:400 v/v), followed by five washes in PBT, 
1 wash in 0.5 X PBT, 0.5 X hybridization buffer (1 X hybridization 
buffer is 50% formamide, 5X SSC, 5X Denhardt’s, 0.1% Tween, 
100 pig/ml Salmon sperm DNA), and then 1X hybridization buffer. 
Incubation continued at 75°C for 30 min, followed by prehybridiza¬ 
tion overnight at 60°C. DIG-labeled probes at a concentration of 50 
ng/ml were added to Artemia in fresh 1X hybridization buffer, and 
hybridization was carried out overnight at 60°C. Hybridization was 
followed by (1) three washes in hybridization buffer at room tem¬ 
perature; (2) five washes 45 min each at 60°C in hybridization 
buffer; (3) two washes in 2X SSC, 0.1% SDS for 10 min each at 
60°C; (4) two washes in 2X SSC, 0.1% SDS for 15 min each at 
60°C; and (5) three washes in PBT for 10 min at RT. The animals 
were incubated overnight in anti-DIG antibody (1:2000 in PBT, 2% 
bovine serum albumin), rinsed three times in PBT, washed six times 
in PBT for 10 min each, five times in PBT for 30 min each, and then 
twice for 10 min each in alkaline phosphatase development buffer 
(100 mMNaCl, 50 niM MgCl 2 , 100 mM Tris, pH 9.5, 0.1% Tween). 
Color substrate (5-bromo-4-chloro-3-indolylphosphate [BCIP] and 
nitro blue tetrazolium [NBT]) was added and animals incubated un¬ 
til staining was complete (30 min to several hours). Animals were 
mounted in 70% glycerol and examined and photographed using a 
Zeiss (Oberkochen, Germany) Axiophot microscope. 

RESULTS 

Identification of an Artemia trh gene 

The strategy for isolating and acquiring sequence informa¬ 
tion for a putative Artemia trh (Af-trh) gene involved first 
isolating genomic clones followed by cDNA clone isolation. 
This approach was taken because before isolation, the devel¬ 
opmental distribution of putative Af-trh RNA was unknown. 
Most bHLH-PAS genes have vertebrate and invertebrate 
orthologs. Between orthologs, the most conserved region by 
sequence is the bHLH domain, and they can be over 90% 
identical in amino acid sequence between Drosophila and 
mammals. Therefore the bHLH domain of D-trh was used as 
a probe to screen an Artemia genomic library. 

Two bHLH-positive genomic clones were isolated, and 
sequence analysis indicated the presence of a gene highly re¬ 
lated to D-trh and Bm-trh. The sequence from the Artemia 


gene was used to design PCR primers for RT-PCR. Prelimi¬ 
nary experiments indicated that Af-trh transcripts were 
present at all stages of larval development and in the adult 
(data not shown). PCR of larval Poly(A) + RNA resulted in 
amplification of a cDNA clone of 1.5 kb that included se¬ 
quences encoding the bHLH, PAS-1, and PAS-2 domains. 
Comparison of the Artemia protein sequence with the se¬ 
quences of Drosophila and Bombyx Trh proteins and with 
the human Neuronal PAS1 (NPAS1) (Zhou et al. 1997) and 
Neuronal PAS3 (NPAS3) (Brunskill et al. 1999) Trh-like 
proteins revealed high sequence identity (Fig. 1). The NPASJ 
and NPAS3 bHLH-PAS genes are prominently expressed in 
the CNS. The strong sequence conservation among these five 
proteins includes the bHLH, PAS-1, and PAS-2 domains, as 
well as the spacers between bHLH-PAS 1 and PAS-l-PAS-2. 
The sequence C-terminal to PAS-2 is unconserved, a feature 
commonly found among bHLH-PAS proteins. 

Comparison of the Artemia protein with all other bHLH- 
PAS proteins demonstrates that it is an Artemia Trh (Af-Trh) 
ortholog. Within the bHLH domain, Af-Trh shares 91% se¬ 
quence identity to D-Trh and Bm-Trh but only 70% identity 
to Drosophila Sim (Table 1). Of the Drosophila bHLH-PAS 
proteins that dimerize with Tgo (Sonnenfeld et al. 1997), in¬ 
cluding Sim, Similar, and Spineless, all show significantly 
less sequence identity to Af-Trh than does D-Trh, and com¬ 
parison of Af-Trh with other bHLH-PAS proteins, such as 
human ARNT and Clock, show much lower sequence iden¬ 
tity (Table 1). Similar relationships are also found when 
comparing PAS-1 and PAS-2 domains among bHLH-PAS 
proteins (Table 1). Examination of a phylogenetic tree that 
compares the full array of bHLH-PAS full-length protein se¬ 
quences (Fig. 2) demonstrates that Af-Trh clusters with other 
Trh proteins in a distinct subfamily of bHLH-PAS proteins 
that include D-Trh, Bm-Trh, and the mammalian NPAS1 
and NPAS3 proteins. As expected, the insect D-Trh and Bm- 
Trh proteins are more closely related to each other than to 
Af-Trh. Similarly, D-Trh and Bm-Trh are closer in sequence 
identity to Af-Trh than to mammalian NPAS1 and NPAS3. 

Comparison of the bHLH, PAS-1, and PAS-2 domains 
between Af-Trh and D-Trh show sequence identity of 91%, 
58%, and 56%, respectively (Table 1). The lower conserva¬ 
tion between PAS domains compared with the bHLH do¬ 
main is a relationship also found when comparing other 
bHLH-PAS proteins. The PAS-1 and PAS-2 domain com¬ 
parisons among the Trh subfamily (Fig. 1) differ with respect 
to blocks of sequence identity. PAS-2 (102-103 amino acids 
in length in the proteins compared) has relatively uniform 
distribution of sequence identity scattered throughout the do¬ 
main, with the longest stretch of sequence identity three res¬ 
idues. In contrast, PAS-1 (129-159 amino acids in length) 
begins at an uninterrupted cluster of sequence identity (9 amino 
acids long), followed closely by another stretch of identity (18 
amino acids), and has another block of identity (10 amino ac- 
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Af-Trh 
D-Trh 
Bm-Trh 
H-NPASl 
H-NPAS3 


|--—basic-helix-loop-helix---| 

(1) MRKEKSRDAARSRRGKENYEFYELAKMLPLPPQITSQLDKASIIRLTIAFLKLKEFTSNG 

(77) ..... * »RD*SGH« 

(71) l*»*******************»»»******aa****»****»»****sy«*«rd»sgh* 
( 45 ) q******n***»»****»L»»*»*»»L»*»*GA*SI*******V*»SVTY»R*RR*AAL* 
(21) ...... SY**MRD*ANQ* 


Af-Trh 
D-Trh 
Bra-Trh 
H-NPASl 
H-NPAS3 


|--PAS -1 - - 

(61) DPPWNNMGSGTGK-GLKIHSSLHPP-GFSPDQIEQFQGTHILQSLDGFAFSLGHDGR 

(137) ••••TREA*SSS*-LKSAAIRR-SPAV»LF**H»•••••••••••*LAVAA*•• 

(131) ••••SRDPPPSS*-A* *GAQNRRTAA* LAMELF* »H*••••••••••••L*VAA»•• 

(105) A»•*GLRAA* PPA-••APGRRG--•AALVSEVF**HL*G*••••••••V*A*NQE »K 

(81) •••••LRME* PPPNTSVKVIGAQRRRS•SALAIEVF*AHL* S.*•V*A*NQE *K 


---PAS-1- 

Af-Trh (116) FLYISETVSIYLGLSQVELTGSSVFDYIHTADHQELADQLGINLNLSNSTSSNFPPSPAS 

D-Trh (18 9) ••••••••••••••••••••I••••*Q*•*S•I••••*LS *TSGGGGGGGSSS *GGG 

Bm-Trh (187) . 

H-NPAS1 (159) ••••••••••••••••••••••••PG**S*VLE»*•LRTTTPGPPTPSSVS*SS- 

H-NPAS3 (141) . V*PG**V«M*E***MK*PPGRGLL*QGTAEDG- 

---PAS-1--- 

Af -Trh (176) - GSGN-GDDLGSESTGMSNMSISKDADYKGFDRAFCLRMKSTLTK 

D-Trh (249) GGGAGGGMASPTSGASDD*SGTHGTNNPDVAA*MTQA*T-SG*»*Y**S*»V»*»»***» 

Bm-Trh (246) ---SEE*S-QHGTNNPDVSSQMSLAASGSL*R*M**•••V** •••••• 

H-NPAS1 (218) •-SSSSLADTPEIEA»LTKVPP--SSLVQE*S*FV*****••• 

H-NPAS3 (200) ---A*SASSSSQSETPEPVEST»P*LLTTDNTLE*S*FI* ••• •••• 


Af-Trh 
D-Trh 
Bm-Trh 
H-NPASl 
H-NPAS3 


-PAS-1-] 

(219) RGCHMKSS----GYRVILIVGRMRAQQTGAT S SDG 

(308) • • • • F • • *GYRASDATSNCNNGNNASNNAKNVKNPGSN* S•V•LLCKL•P•Y* FSH* RKS 

(290) ••••F*** ---••••V*MLC*L*P*YSFSH»RKS 

(258) ••L*V*A*-**K» *HVT* •-HA 

(244) • *V • I • • •---••K**H*T**L*LRVSLSHGRTV 


Af-Trh 
D-Trh 
Bm-Trh 
H-NPAS1 
H-NPAS3 


|.- - PAS - 2- 

(250) TKSSKLLGLVALALALPPPAAHELKLESDMFVMRLNFEFKWHCEP-RISDLMEFTSEEL 
(368) Q--PP***M****I****»SV»*IR*»C*»**T*V**DLR»A»»**-*V***LDYSP*D* 
(321) P--TV***M****I*»**»SV»»IR*»*»**»T»I**D*RIA»***S*V*E*LGY*A»»* 

(281) L-•L***GHT»**APLA**P*HGH*I*F**SLGLTILA**S-*V»»H»DLGPS** 

(275) P--•QIM**»W*H*****TIN*VRIDCH***T*V*MDLNIIY**N-••••Y»DL*PVDI 


Af-Trh 
D-Trh 
Bm-Trh 
H-NPAS1 
H-NPAS3 


-PAS-2---- 

(309) TGKSLYSFSHGEDLEKLRRCHCDLINKGQAMSGCLRMLNRHGGFSWVQLCATWCSSKNS 
(425) •LC»A»*ANR**KS*S***E»••VLT*YY»LM»KS**YT*L*T***••••T**A 

(379) •••N»*TLC****AN***K**L»*M**»»VLTHYY*IM»KL**YT*M»T**********A 

(334) V*R*OQ*V**Q*ATRI*QS*V»*LD*»•V»T»YY*W»Q*A***V*L*SV*•*AG*G»SP 
(332) V*»RC*H*I•A**V*GI*HS*L» *L*•••C r VTKYY*WMQKN»*YI•I*SS»•IAINA**A 


Af-Trh 
D-Trh 
Bm-Trh 
H-NPASl 
H-NPAS3 


-PAS-2-| 

(368) EEQSIVCVTYVLSRPDYGGFPMDTSQLETDIKATDIKVETPISETSAPTAEVINICDSGV 
(485) D* *N*I•»N»•I»NRENENMIL*CC*••PSPDSIKHEEGLGNDKS•GS PGGDASGEGNSH 
(439) ••I••N*•I*GRE»PNTV*•CC*•*ESVQGVKREETTGDP*NGP*DTDNTGGPPPTS 
(394) G*HHVLW*SH***QAEG*QT*L*AF»*PASVACEEASSPG*EPTEPE*PT*GKQAVPAEN 
(392) N*KN*IW*N*L*•N»E»KDT»••IA*•PHLPEK*SESS•*SD**SDSKDTSG*TEDNENS 


Af-Trh 
D-Trh 
Bm-Trh 
H-NPASl 
H-NPAS3 


(429) TVEMRPQFNQVQDLDSVEVPTEVARIGSSENEETYAARLDSARESEIVRSPSGVSRRGR 
(545) LSAGDMKL*SPKTDSEGHSHRGRG*SAAASHGSSMNSLTMIKDSPTPLGVEIDSGVLPT 
(499) RHPAERTEPPTNSETASSAPVNASDVTHLTRLNDAQPLPLHTTPERTSPM*ARRLKKRK 
(454) EAPQTQGKRIKVEPGPRETKGSEDSGDEDPSSHPATP»PEFTSVIRAGVLKQDPV*PWG 
(452) KSDEKGNQSENSEDPEPDRKKSGNACDNDM»CNDDGHSSSNPDSRDSDD*FEHSDFENP 


Fig. 1. Protein sequence of the Artemia Trh protein is highly related to other Trh and Trh-like proteins. The Af-Trh sequence is shown on 
the top line and spans 487 amino acids. It likely has additional residues both N-terminal and C-terminal. The Af-Trh sequence is aligned 
to the Trh proteins of Drosophila (D-Trh; NP_523872) and Bombyx (Bm-Trh; BAA22946) and the Trh-like human Neuronal PAS1 
(H-NPASl; AAB47248) and Neuronal PAS3 (H-NPAS3; NP_071406) proteins. The alignment was determined using CLUSTAL W. The 
sequences after Af-Trh residue 395 are shown without gaps because identity is low. •, Sequence identity to Af-Trh; -, Gaps. Each se¬ 
quence is numbered on the left. At the top are shown the extents of the basic-helix-loop-helix, PAS-1, and PAS-2 domains. 
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Table 1. The percent sequence identity between Af-Trh 
and selected insect and mammalian bHLH-PAS 
proteins is shown quantitatively for the 
bHLH, PAS-1, and PAS-2 domains 

Gene 

% Identity to Af-Trh 

bHLH 

PAS-1 

pas-: 

Bombyx Trh 

91 

55 

60 

Drosophila Trh 

91 

58 

56 

Human NPAS3 

85 

50 

42 

Human NPAS1 

72 

48 

44 

Drosophila Sim 

70 

29 

37 

Human HIF-la 

66 

24 

43 

Drosophila Sima 

62 

25 

30 

Drosophila Ss 

37 

22 

22 

Human Ahr 

35 

23 

19 

Human ARNT 

26 

18 

23 

Human Clock 

24 

18 

29 


Sequences were aligned using ClustalW. The highest % identity 
was observed in all three domains between Af-Trh and the Trh pro¬ 
teins from the insects. Bombyx and Drosophila, indicating that the 
Artemia gene encodes a Trh homolog. Sim, Single-minded; HIF-la, 
Hypoxia inducible factor-la; Sima, Similar; Ss, Spineless; Ahr, Aryl 
hydrocarbon receptor; ARNT, Aryl hydrocarbon nuclear receptor 
translocator. 


ids) near the end. Also contained in PAS-1 are two large 
blocks of low sequence identity and two insertions in D-Trh 
of 15 and 29 amino acids not found in the other sequences. 
These differences may reflect PAS-1 and PAS-2 species- 
specific biochemical interactions and functions and different 
functions between PAS-1 and PAS-2. 

Af-trh is expressed in the larval salt gland 

Expression of the Af-trh gene was examined in developing 
larvae by in situ hybridization (Fig. 3). The early naupliar 
stages showed high levels of trh expression specifically in 
the salt gland (Fig. 3, A and C). The osmoregulatory salt 
gland is one of the most prominent morphological features of 
the nauplius, because of its prominent position on the dorsal 
surface of the head region. Salt gland expression was ob¬ 
served in emerging stage I naupliar larvae (Fig. 3A), indicat¬ 
ing that Af-trh expression is present in the pre-nauplius when 
the salt gland is forming. Expression continued in the salt 
gland during further naupliar development (Fig. 3C). The 
early expression pattern of Af-trh mRNA is consistent with 
it playing a role in the development of the salt gland. 

The specificity and validity of Af-trh in situ hybridization 
was confirmed by hybridization with sense and antisense 
probes generated from Af-trh and Artemia engrailed ( Af-en ) 
cDNA clones. Expression of Af-en in segmentally repeated 
stripes of Artemia larvae has been previously reported using 
in situ hybridization (Manzanares et al. 1993) and immuno- 
staining (Manzanares et al. 1996) with the 4F11 En mono¬ 
clonal antibody (Patel et al. 1989). We repeated these results, 


showing that the Af-en antisense probe detects expression in 
stripes (Fig. 3E) similar to immunostaining with Mab 4F11 
(Fig. 3F). More importantly, neither th t Af-en antisense (Fig. 
3E), Af-trh sense (Fig. 3, B and D), nor Af-en sense (not 
shown) probes showed expression in the larval salt gland. 
Furthermore, the sense probes had low background and no 
discernible staining pattern. These control experiments pro¬ 
vide support that Af-trh is expressed in the salt gland. 

Af-trh is expressed in the thoracic epipods 

As the salt gland disappears during larval development, its 
function in osmoregulation is replaced by the thoracic epi¬ 
pods. Correspondingly, Af-trh expression begins to decline 
in the salt gland of the late nauplius and appears in the epi¬ 
pods of the thoracic appendages (Fig. 4, A and C). Negative 
control hybridizations with an Af-trh sense riboprobe do not 
show epipodal staining (Fig. 4, B and D). Artemia thoracic 
appendages develop in a graded manner with the anterior 
segments developing before the posterior ones. Figure 4A 
shows a gradient of Af-trh expression, decreasing from ante¬ 
rior to posterior. The first two thoracic segments strongly ex¬ 
press Af-trh, the third segment has weak expression, and the 
remaining posterior segments lack Af-trh expression. Older 
larvae show Af-trh expression in the epipods of all segments. 
The individual lobes of dissected limbs are shown in Figure 
4C. Localization of Af-trh RNA in these dissected limbs is 
weakly detected in the pre-epipod and robustly detected in 
the epipod. Later in development the Af-trh message is ex¬ 
pressed exclusively in the epipods. 


DISCUSSION 

Using a D-trh bFILH cDNA probe, we isolated an ortholo- 
gous trh gene from Artemia. The conceptual Af-Trh protein 
sequence is 91% identical to the D-Trh and Bm-Trh bHLH 
domains, a degree of identity considerably higher than to any 
other bHLH-PAS protein. Similarly, the PAS domains of Af- 
Trh are significantly more related to insect Trh than to any 
other bHLH-PAS proteins. Analysis of genomic and cDNA 
clones provided evidence for only a single trh- like gene in 
Artemia, although the possibility of multiple Af-trh genes 
has not been rigorously eliminated. The high degree of se¬ 
quence conservation between Af-Trh and insect Trh is not 
surprising, given the high sequence conservation found be¬ 
tween highly related bHLH-PAS proteins of C. elegans. 
Drosophila, and vertebrates. Caenorhabditis elegans has 
several bHLH-PAS proteins highly related to insect and ver¬ 
tebrate bHLH-PAS genes (Consortium 1998; Powell-Coff- 
man et al. 1998; Epstein et al. 2001; Jiang et al. 2001), but 
there is no clear Trh ortholog. However, the C-Pas-A protein 
is in the same phylogenetic branch as Trh and Sim. Its func¬ 
tion and expression pattern have not been reported. Are there 
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Fig. 2. Phylogenetic tree between Af-Trh and other bHLH-PAS proteins. Full-length sequences of representatives of all published bHLH- 
PAS proteins were compared with each other using CLUSTAL W and placed into a phylogenetic tree based on the neighbor-joining pair 
(NJP) method. The results graphically demonstrate that Af-Trh belongs in a Trh subgroup with Bm-Trh, D-Trh, H-NPAS1, and H-NPAS3. 
Af, Artemiev, Bm, Bombyx\ C, C. elegans; D, Drosophila ; H, human; M, murine. Names and references for each abbreviation can be found 
in Taylor and Zhulin (1999). 


vertebrate orthologs of Trh, and do they participate in tubu¬ 
lar, osmoregulatory, or respiratory functions? The mamma¬ 
lian NPAS1 and NPAS3 proteins are both more related to ar¬ 
thropod Trh proteins than to any other protein. In the adult, 
localization of both proteins is highly restricted to the adult 
CNS and was not detected in lungs or kidney (Zhou et al. 
1997; Brunskill et al. 1999). During embryonic develop¬ 
ment, no expression was observed in the developing respira¬ 
tory system, although NPAS3 is weakly present in cells as¬ 
sociated with the submandibular salivary gland and kidney 
tubules. Although not orthologous to trh, the Hypoxia Induc¬ 
ible Factor (HIF) subfamily of bHLH-PAS proteins plays 
prominent and evolutionarily conserved physiological and 
developmental roles in oxygen homeostasis (Nambu et al. 
1996; Semenza 1998; Crews and Fan 1999; Epstein et al. 
2001; Jiang et al. 2001). 

The Af-trh gene is expressed initially in the salt gland and 
later in the epipods. Although evidence is lacking regarding 
the functional roles of Af-trh in Artemia , its expression dur¬ 
ing epipod development and the dramatic roles of Drosoph¬ 


ila trh in salivary duct and tracheal development suggest that 
Af-trh may play a role in the development and/or function of 
the salt gland and epipods. Both Artemia structures function 
in osmoregulation and express genes encoding proteins in¬ 
volved in salt transport, such as the Na + K + ATPase (Conte 
1984; Holliday et al. 1990; Escalante et al. 1995). Cells of 
both tissues possess extensive membrane folds that serve to 
create a large surface area for transport. These cells also have 
a large number of mitochondria. There are some differences: 

(1) the salt gland consists of a single cell type, whereas the 
epipods consist of two distinct cell types (Martin 1992), and 

(2) the salt gland possesses yolk platelets and glycogen gran¬ 
ules, which are absent in epipods (Conte 1984). Thus, al¬ 
though the gross morphology and some properties of the salt 
gland and epipods are distinct, they share a sufficient number 
of features to suggest that some of the same genes may be 
regulated by Af-trh in both tissues. In Drosophila, salt regu¬ 
lation in the larva occurs in the dorsal region of the large in¬ 
testine (Murakami and Shiotsuki 2001; Lengyel and Iwaki 
2002). The large intestine was not reported to express D-trh 



350 EVOLUTION & DEVELOPMENT Vol. 4, No 5, September-October 2002 


\ 










\ 



Fig. 3. The Af-lrh gene is expressed in the larval salt gland. (A) Stage I nauplius was hybridized to an Af-trh antisense riboprobe. The larva 
has just emerged from the cyst (*). Expression is observed in the salt gland (arrow). (B) Stage I nauplius hybridized to a negative control 
Af-trh sense riboprobe showing an absence of hybridization in the salt gland (arrow ) or other cells. (C) Dorsal view of a stage II nauplius 
hybridized to an Af-trh antisense probe. Anterior is to the left. Expression of Af-trh is observed only in the salt gland (arrow). (D) Stage 
II nauplius hybridized to an Af-trh sense probe. No hybridization to the salt gland (arrow) or other cell types was detected. (E) Control 
experiment showing in situ hybridization of stage III larva with an Af-en antisense riboprobe. The expression pattern of Af-en transcripts 
in segmentally repeated stripes (arrow) corresponds to the localization of Af-En protein shown by 4F11 immunostaining. Staining is de¬ 
tected in the cytoplasm. This experiment validates the in situ hybridization technique. (F) Immunostaining of a stage III larva with mono¬ 
clonal antibody 4F11 that recognizes Af-En protein. Af-En is localized to nuclei in segmentally repeated stripes at the posterior border 
of each segment (arrow). (A-E) Artemia larvae were subjected to in situ hybridization w ith DIG-labeled riboprobes followed by reaction 
with alkaline phosphatase-labeled anti-DIG antibodies. NBT, and X-phosphate. (F) Larva was stained with En 4F11 antibody and 
treated with horseradish peroxidase-conjugated anti-mouse antibodies followed by diaminobenzidine/H : 0 ; histochemistry. 
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Fig. 4. The Af-trh gene is expressed in epipods. (A) Horizontal view of a stage VII larva hybridized in situ to an Af-trh antisense ribo- 
probe. Hybridization is detected in the epipods (arrow) of the thoracic limbs. (B) Hybridization to a stage VII larva with a negative con¬ 
trol Af-trh sense riboprobe showed no detectable expression. (C) High magnification view of a dissected thoracic appendage hybridized to an 
Af-trh antisense probe. High levels of expression are detected in the cpipod (arrow) and lower levels are observed in the preepipod (ar¬ 
rowhead). (D) Hybridization using a negative control Af-trh sense riboprobe shows no detectable hybridization in epipods (arrow) or pre- 
epipods (arrowhead). 


(Isaac and Andrew 1996), suggesting that the putative role of 
Af-trh in controlling formation of osmoregulatory organs in 
Artemia is not conserved in Drosophila. Because D-trh is ex¬ 
pressed in the salivary duct and trachea, it was proposed that 
trh might function in tubule formation (Isaac and Andrew 
1996). Although one aspect of D-trh is to control tubule mor¬ 
phogenesis of those tissues, it is not required for the forma¬ 
tion of all tubular structures, because other tubules, such as 
Malpighian tubules, do not express D-trh. The D-trh gene is 
also expressed in the embryonic CNS, a nontubular structure 
(Isaac and Andrew 1996). Similarly, the sites of Af-trh ex¬ 
pression, the salt gland and the epipod, are not tubular. 

The differential control of gene expression by Trh in Dro¬ 
sophila and silkworm may be due, in part, to the use of dif¬ 
ferent Trh-interacting coregulatory proteins in different tis¬ 
sue types. The control of Drosophila tracheal gene expression 
by Trh requires the presence of Dfr, but Dfr is not expressed 
in the salivary duct, suggesting that a different coregulatory 
protein(s) is required to interact with Trh to control salivary 
duct gene expression and formation (Anderson et al. 1995). 
The Bombyx ortholog of Dfr, SGF-3/POU-MI, overlaps 
with Bm-Trh expression in the trachea and also overlaps 


with Bm-Trh in the anterior silk gland (Matsunami et al. 
1998), suggesting that Trh:Dfr interactions may contribute to 
tracheal and silk gland expression and development. The Ar¬ 
temia dfr homolog, APH- 7, is expressed in the salt gland, 
suggesting that it interacts with Af-Trh to control salt gland 
gene expression. Thus, TrhiDfr interactions may be con¬ 
served in diverse arthropods, including insects and crusta¬ 
ceans, although the targets of their action may be different, 
reflecting the different tissue types in which they are ex¬ 
pressed. Expression of APH-1 in the epipods was not re¬ 
ported, so that it is unknown whether Af-Trh and APH-1 
could potentially interact to control epipod gene expression. 
Further evidence that POU-III subclass proteins regulate os¬ 
moregulatory processes is found in C. elegans , where the 
CEH-6 Dfr ortholog controls formation of an excretory cell 
involved in osmoregulation (Burglin and Ruvkun 2001). Be¬ 
cause vertebrate POU-III proteins have also been implicated 
in control of osmoregulation, it has been proposed that this 
is a conserved feature of these proteins (Chavez et al. 1999; 
Burglin and Ruvkun 2001), although the single Drosophila 
POU-III protein, Dfr, is not expressed in the large intestine 
(Anderson et al. 1995), the site of larval osmoregulation. It 
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will be of interest to see whether nematode and mammalian 
POU-III proteins interact with bHLH-PAS proteins to carry 
out these functions. 

The trh gene is prominently expressed in the respiratory 
organs of two insects. Drosophila and Bombyx, and has been 
shown in Drosophila to play a fundamental role in tracheal 
formation. This raises two related issues regarding Artemia 
Af-trh expression in the epipod. Is the function of D-trh in 
development of the respiratory system conserved in Artemia, 
and does the expression of Af-trh in the epipod provide in¬ 
sight into the evolution of trachea? The first issue concerns 
the role of the epipod in Artemia respiration. Respiration in 
small aquatic arthropods typically occurs by diffusion over 
much of the body surface of the entire organism, and this 
also seems likely for Artemia (Martin 1992). In one study, 
chemical ablation of the Artemia epipods resulted in loss of 
osmoregulation, yet a fraction of the treated animals sur¬ 
vived for extended periods of time (Croghan 1958a). This 
has been interpreted as evidence supporting a primary role of 
the epipods in osmoregulation, but not respiration. Although 
the epipod is unlikely to be the principal organ for respiration 
(Martin 1992), it may still contribute to respiration (McMa¬ 
hon 2001), and expression of trh in the epipods as well as the 
trachea offers an intriguing link between these two struc¬ 
tures. In addition, some Crustacea, such as amphipods (Schmitz 
1992), do possess epipodite gills, indicating a relationship 
between epipods and respiration. There are two interpreta¬ 
tions for the occurrence of trh in the insect trachea and bran- 
chiopod epipod. First, the occurrence of D-trh in the trachea 
and Af-trh in the epipod is not evolutionary related but rep¬ 
resents co-option of the gene regulatory cassette for use in 
two distinct tissues. Second, the trh expression in the two tis¬ 
sues is evolutionary related, and a trh -expressing epipod- 
like ancestral structure gave rise to both the insect trachea 
and Artemia epipod. Because of morphological association 
in the fossil record, it has been proposed that the tracheal spi¬ 
racle and the insect wing have co-evolved (Kukalova-Peck 
1978). Molecular and morphological data have supported a 
theory that wings evolved from the epipod of an ancestral 
multibranched appendage that functioned as gills (Wiggles- 
worth 1973; Kukalova-Peck 1983; Averof and Cohen 1997). 
This structural association of the wings and the spiracles to¬ 
gether with the proposed epipodal origin of wings supports a 
model that the trachea and the wings are derived from an an¬ 
cestral epipod. Expression of trh in the epipods and the tra¬ 
chea is consistent with an evolutionary link between these 
two tissues. Identification of trh in other arthropods and 
analysis of its expression and function will provide addi¬ 
tional insight into this issue. 
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